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ABSTRACT 


The Southern California gravity monitoring project begun In 
May 1974, Is Intended to coordinate gravity measurements with the 
long-basellnc three-dimensional geodetic measurements of the ARIES 
(Astronomical Radio Interferometric Earth Surveying) project which uses 
radio interferometry with extra-galactic radio sources. Gravity data 
from 28 of the stations, monitored on an approximately one- to two- 
month basis, have a single-reading standard deviation of 11 Ugal which 
gives a relative single determination between stations a standard 
deviation of lb Ugal. The averaging of data reduces the uncertainty 
and. If gravity does not change during the averaging time. It appears 
that gravity at a station relative to the base can be determined with 
a standard error of 2 to 3 ugal. Where stations could not be placed 
on low -porosity bedrock, the effects of variable groundwater levels 
must be considered. The largest gravity variation observed, 80 Ugal, 
correlates with nearb> water-well variations and with smoothed rain' all 
Smoothed rainfall data appear to be a good Indicator of the qualitative 
response of gravity to changing groundwater levels at other supra- 
sedlment stations, but frequent measurement of gravity at a stai.'on Is 
essential until the quantitative calibration of the station's respo.se 
to groundwater variations Is accomplished. The largest earthquake to 
occur during the survey time near the gravity network was the August 13 
1978 Santa Barbara Channel event (Mj^ » 5. 1-5. 7, ■ 5.6). The closest 

gravity station to this earthquake, 67 km east of the epicenter, also 
exhibits the network's largest gravity change that cannot be related to 
factors other than tectonic distortion. This change is a 50 Ugal low 
occurring from mld-1975 to mld-1977. 
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INT RODUCTIO N 

Changes in the acceleration of gravity at the surface of the 
earth are due to relative movement of mass within the earth, and to 
variation of the separation of the measurement point from the earth's 
center of mass. Atmospheric effects (discussed by Warburton and 
Goodklnd, 1977) are significantly smaller than other uncertainties in 
local gravity surveys and, for our purposes, are not con- 
sidered. Because the processes of tectonophysi cs can alter both 
density of subsurface rocks and the elevation of the ground surface, 
changes in gravity can be an important geophysical tool for understanding 
the forces and deformation of the earth's crust. 

Gravity measurements are essential complements to measurements of 
elevation changes In order to Interpert the sub-surface crustal 
distortion taking place In tectonically active regions. Elevation 
measurement programs now In progress measure either the geometric 
elevation change directly, as In the extraterrestrial methods such as 
Very Long Baseline Inferometry (VLBI), or they measure the orthometric 
elevation referenced co the geold which Itself may change (for a 
discussion, see Whitcomb, 1976). The program described here was begun In 
order to coordinate gravity measurements with VLBI project ARIES (Astronomical 
Radio Interferometric Earth Surveying) whose survey sites have spaclngs 
of 100-1000 km (see, for example, MacDoran, 1974). The gravity network 
was also densifled to Intermediate sites to provide better spatial 
resolution between the ARIES survey sites and to correlate with leveling 
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dat.i In tiu* Mnini' rogloMH. 

Om* of tht* ill ff lout t li‘H with most geodetic technlqucH 1 h that 
their expeiiMo precltulen high ri‘Holiitlon In both apa^lal and tem|>oral 
coverage of a region that Is suspected of rapid tectonic distortion, 
la'vellng surveys give detailed coverage along the line of survey hut 
are generally repeated only every few years. Kxt rater rest rl al 
techniques such as VhHI or laser ranging currently are n»>st useful 
at station distances ol 100 km or more and, with a few exceptions, are 
repeated vearlv at best. Gravity meter surveys are relatively In- 
expensive compared to the above methods and therefore provide .1 
means for bridging the temporal and spatial gaps left by widely- 
spaced or Infrequent surveys to better define the time and extent ol 
a crustal distortion event. To best corapleim’iU other geodetic 
surveys, a gravity survey must therefore be both widespread to cover 
the spatial spread of extraterrestrial techniques and be frequent to 
cover the temporal gap left by surveys that are too expensive to 
repeat at short time Intervals. These considerations have dictated 
the nK>de of the gravity survey ilescrlbed here. It was designed to cover 
a wide geographic range, to repeat every one to two months, to co- 
locate with as nunv geodetic and geophysical me.isuremt'nt sites as 
possible, and be cost effective. 

Tlu' first consideration In the developmi*nt of a time-dependent 
gravity monitoring program In tectonophyslcs Is to determine the 
expected changes In gravity that are related to crustal distortion. 

Tlie so-called "I’almilale Bulge" In Southern California (Castle et al., 

147b) represents a broad elevation ch.ange of up to em over approxlmat«*lv 
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10-15 yearti. Murh of thla occurred during tiroes at least as short 

an the interval between leveling surveys, of the order of 3 to 5 years. 

More recent evidence indicates elevation decreases of 17 cm over a period 

of two years or less north of Los Angeles. If all of this elevation 

change were to take place with no horizontal mass transfer, then the 

imwement will be characterized by the free-air distortion gravity gradient 

of 3.08 Ugal/cm. In this case, a 17 cm elevation decrease would give 

a gravity increase of 52 hgal and a 45 cm increase would give a 

gravity decrease of 139 pgal. Density changes within the earth's 

crust can also strongly affect gravity. Gravity and elevation data 

from the Mataushlro earthquake swarm in Japan (Klsslinger, 1975) 

show gravity changes of up to 80 hgals wliich liad up to 40 Ugals 

deviation from tlie free-air distortion gravity gradient. Tills 

-4 

combined gravity and leveling data indicated a 0.6-1. 8 x 10 
dilntant volumetric strain within the crust (Whltcomli, 1976). During 
a 3 to 4 month period prior to the M ■ 7.8 Tangshan, China, earthquake, 
gravity changes were measured ranging from -90 to +167 hgals with 
the positive change occurring over the region of the future shock 
(Wang, 1979). Because no coincident elevation data were available, 
this 257 Vigal range could have been the result of either elevation 
change or subsurface density change. For example, gravitational effects 
of this size are well within the estim.ited range of density variation, 
due to dllatancy models of a crust in preparation for a major earth- 
quake (Wliltcomb, 1976). 

As will be shown here, the estimated standard deviation of 
.1 meter reading at a station Is about 11 hgal . Since single gravity 
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detormlnat lonH In the* survey are relative to a reference and are not 
absolute readings of gravity, each determination must consist of at 
least two readings, one for the station and one for the hase. Thus, 
a single determination should have a combined standard deviation of 
16 Ugal. Tlie averaging of repent determinations can reduce the 
standard deviation to some extent provided that there are no systematic 
errors and that gravity does not change between readings. Two standard 
errors of the mean of a five-point moving average used here is 
calculated to be about 10 bgal, and this value is considered to be a 
reasonable estimate for testing the significance of a change of the 
nK>vlng average. If 10 Ugal is translated into a free-air distortion 
elevation change, it would be equivalent to 3.3 cm of elevation change. 

If the 10 pgal is translated into a pure subsurface density change 

for a body similar to the size of the Matsushlro distortion, thickness 

of 5 km and radius of 5 km, it would be equivalent to about 3 x 10 ^ 

dllatant volumetric strain (Wliltcomb, 1976). Tlie est invited 

accuracies of the gravity survey are indeed adequate to provide significant 

insight into tect ononhvs I cal processes. 

The purpose of this paper is to describe the Southern 
California gravity survey and data reduction, show the temporal 
variation ot gravity data since 1974, and show possible correlations 
with rainfall, well levels, and earthquakes in the region. 
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GRAVITY MEASUREMENTS 

The gravity Rtatlon network as listed In Tahle I and shown in 
Figure 1 has Its highest density In the selsnilcally active San 
Gabriel section of the Southern California Transverse Ranges and 
adjoining San Andreas Fault, with outlying stations at the radio 
telescope stations at (kildstone In the Mojave IVsert and (\rens 
Valley to the north. Both of the radio telescope stations are 
used as base stations for the ARIFS long-baseline geodesy program. 

The base station of the gravity network Is circled In the figure. 

Ix)cated at the California Institute of Technology In Pasadena, 

this station Is the starting the ending point for al! surveys and 
all gravity values are direct Iv referenced to this single point. Tills 

Is the sant* base station used by Oliver et al. (1975) in their study of 

the 1971 San Fernando earthquake. 

The gravity meters used for the survey are La Coste-Romherg 

Model G meters, G395 and G465, with electronic readout. 

The earlier part of the survey, from 1974 through 1976, made use 
of G141 which Is the same type of meter without electronic readout. 

The model G Is generally considered to be the state-of-the-art for 
portable gravity Instrumentation with the possible exception of 
the La Coste-Romberg Model D. The Model D has Increased sensitivity 
of the internal screw, one milllgal per revolution versus 70 milllgals 
for the Model G, at the sacrifice of dynamic range. The range of the 
Model D Is 200 milllgals versus 7000 milllgals for the Model G, which 
Is good for the entire earth's surface. In a slde-by-slde survey under 
the same conditions, the Model D appears to be superior only If the 
transport distance Is short and the gravity range Is less than 50 milllgals 
(H. W. Oliver and S. L. Robbins, personal communication, 1975). Survey 
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errors lor the Model G meters used here appear to be dominated by non- 
linear changes in spring length or tares, and this is a characteristic 
that should be shared by Model D. These tares become more of a problem 
with incteased severity of vibrations to the meter due to transport. 

For the type of survey used here, with substantial transport distances 
and ranges of gravity up to 440 milligals, the Model D does not appear 
to be advantageous to the Model G. 

The method of survey is to begin at a base station, follow a 
broad loop of stations, and return to the base station within eight 
to ten iiours for a measurement of closure. Because the long-term 
linear drift during this time is 0 to 13 Ugals, depending on the meter 
used, linear drift should be a minor contribution to the uncertainty 
of the measurement. This method of measurement allows a broad geographical 
region to be tied directly to the value of gravity at the base station 
without accumulation of systematic errors as in the extension of a 
survey line by the Joining of segments. The major benefit of this method 
is that it is the most cost-effective means to monitor gravity over a 
b.oad region. This factor allows more frequent measurements that help 
to establish the time of a rapid change of gravity, and provide a means 
to Improve the accuracy through a combination of frequent repeated 
measurements. 

Closure is the difference In reduced gravity at the base station 
between the measurements at the start and the end of a single survey. 
Statistics of individual meter closures show that two-thirds of the 
closures are less than the following values: G141 - 30 pgals; 

G393 - 25 Ugals; 0465 - 45 ugals. The high value for G465 is due to 
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its young age and presumably will improve after its initial settling 
down period. Overall statistics show that 75Z of the individual 
surveys have closures less than AO gal. On the basis of this number 
and an observation of the general scatter of data with 

higher closures, 40 closure is set as the limit of data acceptance. 

All data associated with higher closures is rejected. 

Very large mlsclosure sometimes occurs and this often has been 
traced to a short episode, of duration less than travel time between 
stations, in which a tare has occurred. The behavlr r of the meter, as 
deduced from the stations' gravity values from previous and following 
surveys, appears to be normal immediately before and after the tare. In one 
such case during a survey in October 1978, G395 was accidentally struck 
which produced a tare of b37 Ugal. Gravity values from the following 
stations in that day's survey showed no unusual deviation from the expected 
values except for the constant 637 Mgal shift at the time of the 
accident. However, meter drift was abnormally accelerated for one month. 

It has been noted in this data set that when large misclosures occur, 
the scatter of individual station gravity values is often larger as 
comapared to moving averages of the data, even when the data reduction 
process assumes a single tare during the day's survey. This may suggest 
that tares tend to occur in groups of two or more, possibly as a result 
of unfavorable meter transport during that particular survey or the nature 
of non-linear drift characteristics of the meter's spring. 

Distribution of the closure error is complicated by the fact that 
closure Is sometimes larger and of opposite sign than what would be 
expected from the long-term linear drift of the gravity meters. Figure 2 
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shows Che drift of G141, 095, and G465 at the base station. The 
amplitude is in dial turns which are approximately 1 milllgal. The 
large shift in G465 in 1978 corresponds to a readjustment of the meter. 
Maximum linear meter drift is about 40 pgals per day and this would 
mean a 17 pgal mlsclosure over a 10 hour survey. Drifts for G141 
and recent data from C )95 arc considerably below this at 4.6 Ugal/day 
and -8 Ugal/day, respectively. Another potential source of mlsclosure 
is from tares or non-linear changes in spring length over short time 
periods. In order to account for both linear drift and tares, the 
mlsclosure distribution is calculated in two ways. First, a linear 
drift with ,ime is assumed; second, a combination of a single tare 
plus a linear drift with time Is assumed. The distribution is 
chosen that produces the best fit, in a least squares sense, to the 
average values of the individual stations. In this manner, the added 
degree of freedom, a Care between stations, is allowed if the gravity 
values for Chat run are closer to the average station values. Thus, 
the extra degree of freedom in distribution of uiisclosure tends to 
minimize any deviation from the average values of Che stations. This 
algorithm chose the simple linear distribution 25X of the time and the 
tare plus linear distribution for the remaining 75Z. 

Errors in the gravity data can be introduced by inaccurate 
estimation and removal of the tides, both solid-earth and ocean tidal 
components. While the solid-earth gravity tides are well knovm for the 
levels of accuracy needed here, the ocean-loading tides can be 
a problem in a region near the coast such as Southern California. In 
some cases errors up to 16 Ugal are possible (Wliltcomb, 1979). The 
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gravity data here have had selected ocean-loading tidal cooponents rsmoved 
as will be described in a separate paper. It is estimated from comparison 
of our solid-earth and ocean-loading tidal calculations to observed tidal 
data that tidal uncertainties are less than +1 ugal with somewhat higher 
errors within 50 km of the coastline. 

In order to compare the results from one meter to another, the 
calibration factor of each meter must be defined more precisely than 
the values provided in the manufacturer's tables, which are at 100 mllllKal 
intervals and given to the nearest 10 Ugal. After the tables are used, 
calibration can be done from the data themselves by computing the statistics 
for each station when multiple meters are used on the same survey. In this 
manner a separate calibration factor is computed for each station that minimizes 
the difference between the meters when they are used simultaneously on surveys. 
Only data with closures of 40 ugal or less are used to compute the calibration 
factors between meters. Because the meter G395 has the longest 
history and is the only meter with data coincident with both 0141 and 
0465, it is used as the reference meter. The individual :itation 
results showing the difference between 0395 and the other 
meters for each station are shown in Figure 3. Data points with no 
error bars, which are standard deviations, are computed with fewer 
than 4 observations for that particular station. All gravity data here 
is adjusted to be equivalent to readings on 0395. 

A potential source of error with a cyclical nature has been 
attributed to non-linearities of the internal screw of model 0 meters. 

The error has a period of approximately 70 mllligals or exactly one 
revolution of the Internal screw. For some meters, this error can be 
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as much as 30 ugals peak-co-peak (R. Jachcns, personal cosmiunlcat Ion 
1978). At Che currenC drifc race of C141, a complece 70 mllllgal 
cycle would Cake abouc 41 years. THe lasC cvo years' drifc race 
on G395 would produce a 70 mllllgal cycle In abouC 20 years. The 
polncs In Figure 3 exhlbic characCerlsClcs ChaC mlghc be chls 
Cype of error buc Che daca are noc of a qualify 
as yec Co confirm chls. When Che phase and ampllcude of Che 70 
mllllgal-perlod varlaclon is calibraCed, Chen Che daCa can be correcCed 
If necessary. FurCher sCudy of Chls behavior la In progr<.js. 

Figure 4 shows Che gravlCy daCa as a funcClon of Clme for 28 
sCaClons of Che Souchem California area as shown In Figure 1. All 
gravlCy daCa here represenC Che gravlCy aC Che scaclon minus Che 
gravlCy aC Che Pasadena base sCaClon. DlfferenC symbols refer Co 
Che mecer used for Che reading; Crlangles for G141, squares for G395 
and circles for G46S. The scaclon name and number along with Ics 
average gravity value are shown In each diagram. The average value 
has been subtracted from the data for each station so that the zero 
line represents Che average. The error bars on Individual points In 
Figure 4 are estimated solely from the closure. The total length 
of the error bar la the closure or 20 Ugal, whichever the larger. 

Thus, because the closure limit Is 40 Ugal, the error bars for a 
single point range from ±10 to ±20 Vigal* Assuming ChaC smaller 
mlsclosure Implies a lower uncertainty for Che data points associated 
with that survey, these error estimates are in good agreement with the 
16 ugal standard deviations estimated below for all data. 

The solid lines In Figure 4 represent ±2 standard errors of the 
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mean of a moving, flva-polnc weighted average of the lndl*d.dual data 
points shown In the saise figure. The average Itself, which Is not 
plotted. Is mldeay between the solid lines. Uelglitlng Is done with 
the assumption that the error bar shown for each point represents that 
point's standard deviation. The i 2 standard error range Is 
used as a test of the significance of any change observed In the 
average. 

Independent confirmation of the reality of gravity changes outside 
the ±2 standard error range Is fortunately available for Station 22. 

Three Independent gravity surveys have been made to this station by the 
National Geodetic Survey. Each survey consists of four to six determinations 
using meters different from those described In this pr.t 'i The data are 
shown as solid points for Stations 22 and 23 In Flgiue 4. The first two 
points show no change while the error range shows a significant Increase 
during mld-1973 and then a return to the early-1975 values by the time 
of the second NGS point In mld-1976. When the moving error range again 
Increased by about 35 pgal In early-1977, the NGS was able to repeat the 
determination and found the same Increase In gravity. Thus the changes 
seen for Station 22 In Figure 4 of the order of 30 pgal are considered 
to be real and take place over periods of 6 months or less. It Is clear 
that In order to avoid temporal aliasing of the data as occurred with the 
first two NGS data, measurements must be made at Intervals significantly 
less than 6 months. 

The scatter or standard deviation of Individual gravity measure- 
ments can be estimated with two essentially Independent methods. The 
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first Is to estimate the standard deviation of data for stations 
close to the base station under the assumption that their gravity 
relative to the base Is not changing with time. At least two stations 
fulfill these criteria. Stations 16 and 24. The standard deviations 
of the data for Station 16 are as follows: G141 - 14 pgal, 

and G395 - IS ugal (G465 does not have enough data). The standard 
delvatlons • f the data for Station 24 are as follows: G141 • 15 Ugal, 

G395 - 14 Ugal, and G465 - 22 Ugal. The readings for these stations 
dre generally not done on the same day and therefore they should be 
Independent data sets. The larger standard deviation of G465 Is presumably 
due to the first year'o Instability of a new neter and a small sample. 
Because two readings are required for each point, the base station and 
the survey station, single-reading standard deviations for G141 and G393 
would be 11 Ugals and for G465, 16 Ugals. A second method of estimation 
of a single-reading standard deviation Is to use the comparlsun of two 
metjr readings for the lame station and survey, which should be Identical 
after the meter calibration factor Is removed. The standa'i.'* deviation 
of 85 observations of G141 compared to G39S Is 22 ugal. The standard 
deviation of 79 observations of G465 compared to G395 Is 23 ugal. Each 
of these comparisons requires four single readings, one base station and one 
survey station tui <;ach meter. Thus, a single-reading standard deviation 
for each comparison is about il ugals. In agreement with the independent 


estimate above. 
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!«• s p ON s F To^ * 

Although every attempt haa been made to place gravity stations 
at bedrock sites where near-surface porosity Is low. It Is not always 
possible to avoid sediment basins. This Is due to factors such as 
the absence of bedrock at a tectonically Important site, ease of 
access, and the fact that many survey bench marks are located In valleys 
along maior roads. One half of the sites are situated on sedlmi‘nts. 

The gravitational effect of a varying level of groundwater can be significant 
but It depends on factors such as porosity and extent of the aquafer. The 
attraction due to an Infinite horizontal slab of water Is 0.419 Mgals 
per centimeter of slab thickness. The general behavior of groundwater 
levels depends on the general rainfall of the region with modification 
by nearby pumping If present. The presence of pumping uswilly ensures 
a program of groundwater level monitoring that can be used to 
directly estimate the gravitational effect of level variations. However, 
because mi')nltorlng wells are usually not located lirant'dlately next to a 
.gravity station. It Is found that well levels do not always correspond 
to gravity variations at a station. For the station set, a search 
was made for monitor-well data c.i-e to all gravity stations located 
on sediments. Tlie gravity data for these stations can then be correlated 
with their respective well data. In addition, data for all stations can 
be correlated with precipitation, tfte ultimate source of aquafer 


recharge. 
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Figure 5 shows monthly rainfall In Los Angeles (Station 716) 
after it has been passed through a lo%r>pass exponential filter with 
a drop-off of 1/e in one year (the filter was started in 1959). 

Treated in this manner, the rainfall data is intended to resemble the 
behavior of an aquafer that is chr.gt>d with water during the winter 
rain season in Southern Calflomia and slowly drained either by lateral 
flow or by pumping. The average rainfall is 1.27 Inches/month, but 
statistically, 76Z of the rain comes ir the months of December, January, 
February and March. The figure shows that rainfall in 1974-75 was 
average, 1976-77 was low, and 1978-79 was high. 

The largest gravity change of the stations in Figure 4 is that 
of Station 28, which shows a 80 ugal increase in early 1978. If this 
is compared with the filtered rainfall data of Figure 5, a clear 
correlation 1s seen. During the heavy rains of early 1978, the closest 
monitored well to Station 28, Well 7128C, 100 meters away, showed 

a very large increase in level of nearly 16 meters. If this well data 
were to be converted to a gravity change with the assumption of a 
horizontal water table in an aquafer of porosity 0.13, it would closely 
fit the gravity data of Station 28 as shown in Figure 6. 

The zero level in Figure 6 is arbitrary. Thus, both the amplitude 

and the timing of the large gravity increase at Station 28 is readily explained 

by groundwater. However, for this station, the converted groundwater data of 

Figure 6 in 1975 are significantly higher than the measured gravity 
of Figure 4. This may be due to a different phase response 
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of the aquafer under the gravity station compared to the well. For 
example. If the smoothed rainfall data of Figure 6 Is compared to the 
gravity data, It can be seen that a slightly faster dropoff of the 
early 1975 rainfall charge could bring the late 1975 levels down to 
the 1976 levels, as Is seen In the observed gravity. Tlius, in the case 
of the highest gravity response to groundwater variations, smoothed 
rainfall appears to be a good Indicator of the qualitative response 
of a gravity station to changing groundwater levels. 

Figure 6 also shows the converted groundwater data for all other wells 
that are near sedlment-basln stations and for which water level 
monitoring data has been found. The well number and corresponding gravity 
station name and number are shown for each data set. The well level changes 
were all converted to gravity change with the same assumptions of a 
horizontal groundwater slab In an aquafer of 0.13 porosity. Data 
from well 4076 Is close to the base station (CALTECHL) and shows changes 
that would correspond to a drop of 20 pgals during 1976-78. Two 
closeby gravity stations, 16 and 24, are In sites that should have 
little groundwater effect and were Installed to mi>nltor the behavior 
of the base station. Gravity data for Stations 16 and 24 In Figure 4 
show no Increase In gravity during the 1976-78 period during the low 
of data from well 4076 or the low of the filtered rainfall. It is 
concluded that the base station has had no measurable response to 
nearby well data or rainfall as yet. Well 5873D close to Station 27 
shows little change, yet the gravity data have a good correlation with 
the smoothed rainfall data of Figure 5. Well 8488A close to Station 29 
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shows a large water level variation that correlates fairly well with 
the rainfall data but the gravity data of Figure 4 does not agree. 

In this case the gravity data has a large variation but the low of 
the gravity data Is a few months later than the low of both the rainfall 
and well data. It Is possible that the gravity variation Is due to 
groundwater but that the aquafer system under the gravity station has 
a different phase delay In response to groundwater recharge. Further 
rain cycles should resolve this. The wells 8695A and 8693B near 
Station 31 show little correlation with rainfall and their levels 
have been stable. The Station 31 gravity data has similar behavior 
with no significant change from 1975 through 1978. Well 8876 near 
Station 36 shows little variation In water level but the gravity 
data, which Is scattered after 1978, has some correlation with filtered 
rainfall. 

Some of the gravity stations are on sediments, but no nearby 
monitorlng-well data have been found for comparison. However, the 
above analysis shows that the qualitative behavior of gravltv 
response to groundwater can be estimated from correlation with the 
filtered rainfall data of Figure 5. Of the remaining suprasediment 
gravity stations. Stations 23, 43 and 46 show no rainfall response. 
Station 26 shows some response but the data are scattered. Station 
<»4, located on a bridge benchmark over a stream channel In Wrlghtwood, 
California, shows correlation with the filtered rainfall data but the 
response to the high rainfall In early 1978 is delayed by 3 to 6 
months. Further rainfall correlation study is needed for this station 
to determine Its response. Station 47 shows some correlation with the 
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rainfall but again the response Is delayed. 

Although more data Is needed before an accurate analysis can 
be done, it appears that the gravity stations can be easily calibrated 
as to their response to local groundwater variations by correlation, 
perhaps with phase shifts, with a function constructed from filtered 
rainfall data. The groundwater response can then be removed from the 
gravity data. This analysis points to the need for frequent observation 
of stations in order to compute their gravity response to groundwater 
variations. 

Station AS, on Table Mountain near Wrightwood, California, shows 
a possible correlation with the heavy rainfall of early 1978. Because 
this station is at a high elevation near the top of the mountain, 
heavy show accumulation is possible on the mountain during the cold 
winter and spring months followed by water saturated soils Into late 
spring. Accumulations of precipitation equivalent to 36 cm of water, a 
reasonable estimate for heaw snow pack, that is distributed In a horizontal 
slab under the gravity station would increase the gravity by only 
15 ugals. However, the effect of the downslope water distribution at 
a station on the top of a mountain will increase the vertical gravitational 
acceleration for the samt* thickness of accumulation. A simple 
calculation can be done to estimate the effect. Assume that a 
gravity station Is at the apex of a circular cone as shown In Figure 7 
with sides that slope at an angle ^ from horizontal, _R is the down-slope 
distance, and ^ is the element of cone area at distance R from the 
apex and horizontal angle 0. The mass ^ of a layer with surface 
density p at M is: ^ ■ p dA . The vertical gravitational attraction of 
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Che water layer is therefore: 
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wher'’- G Is the gravitational constant (6.673 x 10 cm g s ). 

Tne first part of equation (1) Is the familiar Infinite slab formula 
so that Che second part In brackets, sin 1 cos 1 In R 2 /Rp Is the 
amplification factor due Co the non-horlzontal mass distribution 
of Che cone. Equation (1) can now be applied Co reasonable values 
for Station 43. If the slope ^ Is 14**, R^ Is SO cm, R 2 Is 1.6 km, 
and the water layer Is 36 cm thick, then the gravitational attraction 
Is 29 Ugals or 1.9 times Che attraction due to a horizontal water 
slab of the same thickness. Half of the contribution of the 29 Ugals 
comes from the first 30 meters of the down-slope distance R. Thus, 
relatively local distributions of snowpack downslope from a hilltop 
site can easily double Che attraction compared to an equivalent 
amount of snovpack distributed horizontally. The value of 29 Ugals 
compares well with the observed increased gravity in early 1978 
at Station 45. One additional mountain site. Station 37, appears to 
have a similar effect although the site Is on a mountain saddle 
Instead of a peak, which reduces the amplification effect. 

The sites along the Malubl coast. Stations 4-9, show no correlation 
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with rainfall with the possible exception of Station 9 at the western- 
most end along the coast. Station 9 shows a 50 pgal low beginning 
In late 1975 and ending In late 1977, a two year duration. This 
gravity low approx! m.itely coincides withe the filtered rainfall low 
of Figure 5 except that the gravity low terminates before the rainfall 
low. This would preclude the possibility of the gravity low being 
directly caused by rainfall Induced groundwater effects. Other 
considerations also weigh against water-related effects. Station 9 
Is situated on a cement structure at the entrance to a highway 
culvert at an elevation above sea level of 3 meters. Because the ocean 
Is nearhy and the ground level Is approximately 2 meters below the 
station on the side opposite the highway, there Is little latitude 
for groundwater- leve 1 variation. Moreover, Station 8 Is in a 
nearly Identical situation and shows no correlation with rainfall, 
therefore, It 1s concluded that the 50 pgal gravity low at 
Station 9 Is not groundwater-induced. 

CORRELATION WIT H TECTONIC ACTIVITY 
Table II and Figure 1 show the earthquakes that have occurred 
In Southern California of Mj -4.5 and larger from January 1974 
through March 1979. Of the events of magnitude 5.0 or larger, only 
four are within 100 km of a gravity station. The first earthquake 
Is a magnitude 5.1-5. 7 (M, ) - 5.6 (M , USGS) that occurred on August 13, 

L« S 

1978 In the Santa Barbara Channel 67 km west-northwest of Station 9. 

As discussed above, gravity data for Station 9 show the only significant 
change among the entire network of stations for which no alternative 
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explanation auch as groundwater variations can be found. Gravity 

decreased by 50 pgals In 1975, Increased again by 1978 after which the 

earthquake occurred as shown In Figure 4 . Under the assumption 

that the gravity change at Station 9 Is caused by tectonic distortion 

and that the tectonic distortion Is related to the earthquake In the 

Santa Barbara Channel, It Is Interesting to compare the station- 

event separation distance with similar precursory distortion data from 

past events. Anderson and Whitcomb (1975) related the maximum 

observation distance fur precursory distortion to the fault 

or aftershock dimension of the earthquake. The relationship closely 

2 

follows the equation (£/L ) ■ -3 where ^ is the length of 

the earthquake fault dimension or aftershock zone and 1. Is the length 
of the affected area of precursory distortion. For the Santa 
Barbara Channel event, the aftershock dimension was approximately 
12 km. This then provides a length L of the distortion area of 
110 km. Thus, the 67 km separation distance of the Santa 'Barbara 
Channel event from Station 9 Is within that which has been observed 
for past precursory distortions. 

If the gravity change observed at Station 9 Is indeed a precursory 
tectonic response to the process of preparation for the August 13, 1978 
Santa Barbara Channel earthquake, then the lack of a similar gravity 
change at Station 8 Imposes stringent constraints on the distribution 
of the distortion area. As mentioned above, the stability of gravity 
at Station 8 provides support to the conclusion that the gravity change 
at Station 9 Is not due to systematic errors for stations in the Malibu 
Coast loop or due to the effects of groundwater. However, because 
the stations are separated bv only 8 km, the distortion area that gave 
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rlHe to Che Rravlty change at Satlon 9 must be either a local effect 
or, if it extends west to the epicenter of the Santa Barbara 
Channel earthquake, its edge must be within a few km of Station 9. 

Either case Is certainly an acceptable physical model, but Che argument 
for association of the gravity change at Station 9 to the Santa 
Barbara Channel event is weakened by the fact that only one station 
clearly shows an anomaly chat can be temporally related to the event. 

The characteristic anomaly duration of earthquake precursors has 
been put into a relation relating the logj^^ T (days) to a linear 
function of magnitude by Tsubokawa (1969), Wliitcomb et al. (1973), 
Myachkln and Zubkov (1973), Scholz et al. (1973), and Kikitake (1975). 

If the larger of the estimates of the magnitude of the Santa Barbara 
Channel earthquake were to be applied to these formulas, they would 
give anomaly duration times ranging from 240 to 420 days. The anomaly 
time of the data in Figure 4 for Station 9, as measured from 
the beginning of the anomaly to the event In August 13, 1978, is 
about 1,000 days. Thus the gravity anomaly duration is longer than 
that expected from the earlier formulas. However, these formulas 
will require considerable testing and probable modification before 
they can be applied confidently, especially for varying types of 
earthquakes and tectonic regions. 

The next two events occurred close to one another In time and space 
on October 4, 1978, 50 km north of Stations 40 and 41. The first was 
a magnitude 5.8 (Mj^) ~ 5.2 (M^* USGS) and the second was 5.3 (M^ ) within 
an hour of the main shock. Tlie time of events is shown by arrows in 
Figure 4, Stations 40 and 41. Unfortunately, due to logistics probl€*ms. 
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only four setn of gravity data were gathered for these stations so 
that temporal details of the behavior of gravity are not available. 

The latest points were obtained Just after the earthquake, and show 
no consistent change from the data taken In 1975 and 1977. However,^ 
it Is entirely possible that a change in gravity could have t.iken 
place during the 413 days between the survey in 1977 and the event in 
1978. 

The last larger event that is close to gravity stations is the 
January 1, 1979 magnitude 5.0 (Mj) earthquake In the Santa Monica 
Bay off Malibu. This event, as seen In Figure 1, is located 
approximately 17 km south-southwest of Station 5 and south-southeast 
of Station b. Although some gravity change Is seen prior to the 
time of the event In Figure 4, especially at Station 5, the changes 
are of short duration and are Just outside the moving standard error bars. 


SUMMARY 

Gravity Is a diagnostic indicator of the process of tcctonophyslcs 
that can alter both densltv of subsurface rocks and the elevation of 
the ground surface. In order to Interpret surface measurements in 
terms of tectonophyslcal distortion within the crust, both gravity 
and elevation measurements are necessary. Tlie program described 
here was begun in order to coordinate gravity measurements with the 
long-baseline three-dimensional geodetic measurements of the ARIF.S 
project which uses radio Interferometry with extra-galactic radio 
sources. Because gravity measurements are relatively economical, 
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they can be used with greater density both spatially and tenporally 
to increase the resolution of a tectonic event. 

From past crustal distortions and actual gravity measurements 
in an active tectonic area, gravity changes of 50 to more than 200 
Ugal are likely. From the data here, the standard deviation of a 
single LaCoste-Romberg model G reading is 11 Ugal which gives a 
relative determination between stations a standard deviation of 
16 Ugal. The moving average of five readings produces an estimate 
of 10 ugal for two standard errors of the mean which is taken to 
be a measure of this survey's resolution of a change in gravity at 
a typical station. The averaging of all data over the full time 
span of this survey results in an average value of gravity at a 
typical station with standard errors of 2 to 3 ugal . From these 
calculations, the accuracy of the gravity survey is deemed to be 
more than adequate for detection of tectonic distortion within 
the crust. 

Although every attempt has been made to place gravity stations 
at bedrock sites in order to avoid the effect of variable groundwater 
levels, it is not always possible to avoid sediment basins. Where 
possible, water-well level data near the suprasedlment gravity 
stations are used to estimate the effect of varying groundwater level. 
This analysis shows that smoothed rainfall data appears to be a good 
indicator of the qualitative response of gravity to changing ground- 
water levels at a station. The quantitative response, which depends 
cn factors such as porosity and the shape of the water table surface, 
can be estimated from the correlation of gravity with smoothed rainfall 

over several rain seasons. Thus, water well data is not necessary 
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for Che calibration of a gravity station's response to groundwater, 
but frequent raeasurenent of gravity at Che station Is essential until 
this calibration Is accomplished. 

It Is found that most of the significant variations of gravity 
at the 28 stations described here are related to groundwater variations. 
Indeed, Che largest gravity variation observed, 80 ugal , occurred 
at a station In a sediment basin and Che gravity change correlates 
both with well level and rainfall data giving an estimate of 
porosity at chat site of 0.13. 

The largest earthquake to occur within or near the gravity 
network was Che August 13, 1978 Santa Barbara Channel event 

• 5.1-5. 7, ■ 5.6). It Is considered significant that the 

closest gravity station to this earthquake. Station 9 on Che Malibu 
coast, also exhibits the network's largest gravity change Chat cannot 
be related to factors other than tectonic distortion. This change 
Is a 50 ugal low occurring from mld-1975 to mid-1977. However, 

If the gravity change Is relat d to the earthquake, which Is 67 km 
to Che west of the SkSCion, Chen the distortion area must be localized 
or Its edge must be within a few km of Station 9. This reasoning 
follows from the lack of a similar anomaly at the next closest station 
which Is 8 km to Che east of Station 9. The separation distance 
of Station 9 from the earthquake, 67 km. Is within past estimates of 
the size of precursory distortion areas. The anomaly duration, 
about 1,000 days, is longer than that calculated from previously published 
relationships relating the characteristic anomaly duration of earthquake 
precursors to the earthquake magnitude, 240 to 420 days. 
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TABLE I. Gravity Stations 


STATION 

NO. 


LAT. 


LONG. 

GRAVir 

(mllllKals) 

CITl 

4 

34. 

1.200 

118. 

30.330 

25.564 

CIT2 

5 

34. 

2.220 

118. 

38.150 

44.815 

CIT3 

6 

34. 

0.100 

118. 

48.350 

53.185 

CIT4 

8 

34. 

3.200 

118. 

57.800 

79.764 

CITS 

9 

34. 

5.030 

119. 

2.300 

96.271 

LA CO BM 

16 

34. 

12.350 

118. 

10.300 

-52.246 

GOLDST 

22 

35. 

25.660 

116. 

53.300 

-133.282 

ECHO 

23 

35. 

25.660 

116. 

53.300 

-118.192 

LVISTA 

24 

34. 

10.579 

118. 

10.461 

-10.271 

TUJUNCA 

25 

34. 

17.369 

118. 

21.957 

1 

00 

SFERND 

26 

34. 

18.869 

118. 

29.118 

-57.604 

NEWHALL 

27 

34. 

24.064 

118. 

32.547 

-48.506 

MINT 

28 

34. 

24.999 

118. 

28.918 

-50.501 

SLEEPY 

29 

34. 

31.251 

118. 

19.200 

-96.838 

RITTER 

30 

34. 

30.668 

il8. 

14.042 

-124.991 

PALMl 

31 

34. 

32.315 

118. 

6.300 

-126.917 

BLOSSOM 

36 

34. 

31.237 

117. 

55.354 

-148.633 

TIESUMIT 

37 

34. 

23.361 

118. 

4.793 

-363.452 

LUCAS 

38 

34. 

17.991 

118. 

9.133 

-154.148 

LACREST 

39 

34. 

13.676 

118. 

11.080 

-88.462 

OVRO 

40 

37. 

13.970 

118. 

16.870 

-116.296 

WESTGARD 

41 

37. 

12.260 

118. 

14.510 

-113.925 

CLYDE R. 

43 

34. 

19.440 

117 

33.828 

-285.611 
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STATION 

NO. 


LAT. 

LONG. 

GRAVITY 
(ml 1 llgals) 

WRIGHTWD 

44 

34. 

21.613 

117. 

38.130 

-343.331 

TABLE MT 

45 

34. 

22.970 

117. 

41.080 

-443.173 

ANTV MT 

46 

34. 

49.110 

118. 

29.820 

-138.691 

ANTV 138 

47 

34. 

46.520 

118. 

36.280 

-147.384 

SKELTON 

48 

34. 

11.800 

118. 

49.300 

-2.011 

OAK FLAT 

49 

34. 

35.957 

118. 

43.187 

-118.901 
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TABLE 

DATE 

197A 

1974 

1975 
1975 
1975 
1975 
1975 
1975 

1975 

1976 

1977 

1978 
1978 

1978 

1979 
1979 
1979 
1979 


II Earthquakos In Southern California of ■ 4.5 and larger . 
January 1974 to March 1979. 




TIME 


LAT. 

U)NC 

■ 


3 

9 

0 

54 

31.91 

34 

23.93 

-118 

28.41 

4.7 

12 

6 

12 

13 

8.31 

32 

42.50 

-115 

23.54 

4.5 

1 

12 

21 

22 

14.84 

32 

45.39 

-117 

59.29 

4.8 

1 

23 

17 

2 

29.43 

32 

57.11 

-115 

29.38 

4.8 

5 

13 

0 

21 

35.58 

34 

59.97 

-119 

6.17 

4.5 

6 

1 

1 

38 

49.23 

34 

30.94 

-116 

29.73 

5.2 

8 

2 

0 

14 

7.73 

33 

j1.19 

-116 

33.48 

4.7 

11 

15 

6 

13 

27.62 

34 

18.22 

-116 

20.48 

4.6 

12 

14 

18 

16 

20.09 

34 

17.38 

-116 

19.30 

4.7 

4 

8 

15 

21 

38.07 

34 

20.81 

-118 

39.34 

4.6 

8 

12 

2 

19 

26.08 

34 

22.78 

-118 

27.52 

4.5 

8 

13 

22 

54 

52.33 

34 

17.31 

-119 

37.58 

5.1 

10 

4 

16 

42 

48.63 

37 

31.68 

-118 

37.89 

5.8 

10 

4 

17 

39 

2.87 

37 

35.06 

-118 

37.04 

5.3 

01 

01 

23 

14 

38.90 

33 

56.70 

-118 

40.90 

5.0 

03 

15 

20 

17 

50.80 

34 

18.30 

-116 

26.30 

4.9 

03 

15 

21 

07 

16.50 

34 

19.50 

-116 

26.60 

5.2 

03 

15 

21 

34 

25.50 

34 

20.80 

-116 

26.90 

4.5 

03 

15 

23 

07 

58.90 

34 

19.60 

-116 

26.30 

4.8 


1979 
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FIGURES 

1. Map of the Southern California area showing gravity stations, 
major earthquakes of > A. 5 or larger from January 1974 to 
March 1979, and major mapped faults in the region. 

2. Drift of the LaCoste-Romberg gravity meters G141, G395, and 
0465 at the base station CALTECHL . The amplitude is in dial 
turns which is approximately one milligal. 

3. Station calibration factors showing the difference between 
G395 and the other meters for each station as a function of 
the station's average gravity value. Data points with no 
error bars, which are standard deviations, are computed with 
fewer than 4 obser «rat ions for that particular station. 

4. Gravity data as a function of time for 28 stations of the 
Southern California area. Stations are shown in Figure 1. 
Triangles refer to meter G141, squares for G395, and circles 
for G465. All gravity values represent the gravity at the 
station minus the gravity at the base station CALTECHL. The 
station name and number along with its average gravity value 
are shown in each diagram. The average value has been sub- 
tracted from each point. The double lines represent plus and 
minus two standard errors of the mean of a five-point weighted 
moving average. Arrows indicate the time of nearby earthquakes. 
Solid points are independently determined values from the National 
Geodetic Survey. 
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5. Monthly rainfall In Loa Angeles (Station 716) after It has 
been passed through a low-pass exponential filter with a drop- 
off of 1/e In one year (the filter was started In 1959). 

6. Water-well level data converted to gravity change under 
assumption of an Infinite slab of groundwater In an aquafer 
with porosity 0.13. The water well numbers are shown along 
with those of the gravity station that Is nearby. 

7. Right-circular cone with slope ^ used to estimate the 
gravitational effect of mountain snovrpack on a mountalntop 
gravity station. 
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